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Abstract 

We consider a model nondispersive nonlinear optical fiber channel with additive white Gaussian noise at large SNR (signal- 
to-noise ratio) in the intermediate power region. Using Feynman path-integral technique we for the first time find the optimal input 
signal distribution maximizing the channel’s per-sample mutual information. The finding of the optimal input signal distribution 
allows us to improve previously known estimates for the channel capacity. The output signal entropy, conditional entropy, and 
per-sample mutual information are calculated for Gaussian, half-Gaussian and modified Gaussian input signal distributions. We 
explicitly show that in the intermediate power regime the per-sample mutual information for the optimal input signal distribution 
is greater than the per-sample mutual information for the Gaussian and half-Gaussian input signal distributions. 


I. Introduction. 

The channel capacity C introduced by Shannon in his seminal work IT] is related to the maximum amount of information 
that can be reliably transmitted over a noisy communication channel. Shannon calculated the capacity of the linear channel 
with additive white Gaussian noise (AWGN) and found the famous logarithmic dependence of the channel’s capacity on the 
signal power: 

C oc log 2 (1 + SNR) , (1) 

where SNR = P/A is the signal-to-noise power ratio, P is the signal power, and A is the noise power. This, in particular, 
means that when the noise power A is fixed, in order to increase the capacity one has to increase the signal power P. 

The interest in nonlinear communication channels has been increasing since the beginning of the 2000’s when fiber optics 
communication systems had to increase both bandwidth and system reach which required the use of ever higher optical power. 
Fiber optic nonlinear channels have been studied both analytically and numerically in numerous papers, see e.g. El, a, a, 
Il5l . 1^ . 171 . l8l . l9l and references therein. The simplified model nondispersive nonlinear optical fiber channel was considered, 
e.g. in ITOl . im . IT 2 I . IT 3 I . im . The investigation of nonlinear communication channels where transmission is affected and 
changed by the signal power is a difficult problem, especially at large SNR 16]. Analysis of the capacity of these channels is 
technically challenging and new techniques and methods are highly desirable to advance these studies El, 03, m, M, ini. 
In this work we consider a simplified model nonlinear channel with a limited range of practical applications. However, methods 
developed for and tested on such model channels might be useful for much more complex and challenging nonlinear fiber 
communication problems. We introduce here a new approach to the calculation of the conditional probability density function 
via the path-integral technique and demonstrate its application using considered model channel as a particular example. 

The channel capacity C is defined as the maximum of the mutual information Ip^ [jfj with respect to the probability density 
function Px[Ar] of the input signal X: 

C= max/p^[x], (2) 

Px [X] ^ ^ J 

where the maximum value of should be found subject to the condition of fixed average signal power: 

P = j VX\X\^Px[X]. (3) 

The mutual information of a memoryless channel is defined in terms of the output signal entropy i2[Y] and conditional entropy 
H[Y\X\. 


Ipx[x]=H[Y]-H[Y\X], 


(4) 
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with 


H[Y\X] = - j VXVYPx[X]P[Y\X]\ogP[Y\X], 
H[Y] = - j VYPout[Y]\og 

Pout 

Pout[Y] = j VXPx[X]P[Y\X], 


(5) 


( 6 ) 


(7) 


where P[F|X] is the conditional probability density function (PDF) for an output signal Y when the input signal is X, 
and Pout\X] is the PDF for an output signal Y. The measure VY is defined as / VYP\Y\X] = 1, and VX is defined as 
f VXPxlX] = 1. The capacity (l2]i, as defined by (l4]i-(|7]i, is measured in units of (log 2)“^ bits per symbol (also known 
as nats per symbol). The input and output signals are functions of time where the signal’s spectrum is restricted to a given 
bandwidth. In general, a sampling of the temporal signal should be introduced to define a discrete-time memoryless channel, 
however, here we consider only per-sample quantities. 

The channel’s mutual information (|4]i depends on the probability distribution Px[X] of the input signal. The input signal 
PDF, that maximizes the channel’s per-sample mutual information is called “capacity-approaching” or “optimal” PDF P^*[X]. 
Obviously, the problem of finding the optimal PDF of the input signal for nonlinear optical channels is of great practical 
importance. 

In the previous studies of nondispersive nonlinear optical channels (e.g. in, m, M) the Gaussian and half-Gaussian 
input signal PDF’s were used as trial functions in order to put low bound constraint on the channel capacity, or to provide 
asymptotic estimate of the capacity in the regime of large SNR. The authors of lfT4l argued, that the half-Gaussian PDF which 
we denote as 


(i)rvi _ exp{-|X|V(2P)} 


Pk^[X] = 


7r|X|(27rP)i/2 


( 8 ) 


provides the best approximation for the “capacity-approaching” input signal distribution at large SNR. In the present paper by 
solving a variational problem we show that it is not the case. We find a true optimal distribution P^*[X] (which in fact is 
different from half-Gaussian distribution) in the regime of large SNR for intermediate power range. We explicitly show, that in 
this regime the mutual information (|4|i for our optimal input signal PDF is larger than the mutual information for the Gaussian 
and half-Gaussian input signal distributions. 

The estimates for the capacity of nonlinear fiber channels with zero dispersion and additive white Gaussian noise in the 
regime of large SNR were obtained in Refs. II3, m. The lower bound for capacity of the channel, based on trial Gaussian 
input signal PDF, reads 113; 


i log (SNR) 


1 


- -/E - log(47r) / log (SNR) \ 
2 ^ SNR ) 


(9) 


where 'yE ~ 0.5772 is the Euler constant. Note that the second term on the right-hand side of Eq. (|9]l was presented as Oil) 
in Ref. lfT3l but it is easily calculated using Eqs. (23) and (24) of Ref. ifTJI . The pre-logarithmic factor 1/2 in Eq. (|9]l arises as 

a result of the fact that in the high power regime, when the signal power P > , the signal-dependent phase noise 

due to self phase modulation occupies the entire phase interval [0, 27r] and, as a result, the phase does not transfer information, 
see Ref. na. Flere 7 is the Kerr nonlinearity coefficient and L is the fiber link length, see below. In ifia capacity estimates 

were also given in the intermediate power range iV <C P ^ . Eor such a power P the following estimate of 

the lower bound for the capacity, based on the half-Gaussian input signal PDE, was derived lfT4l : 


C > — log{jNL) + 


7 ij - 1 -b log(37r) ^ 1 


VySNR 


( 10 ) 


where instead of O the authors presented the explicit function of the parameter SNR which decreases at large SNR, 

see Eq. (40) in M- However, the authors of 03 did not take into account the 1/v/SNR corrections in the output signal 
entropy H[Y], therefore, using these explicit functions in the capacity inequality is beyond the calculation accuracy. It also 
means that the result Eq. (40) of 03 is not a lower bound on the capacity. It is worth noting that in their result there is term 
log 2 missing. Also their result does not recover the Shannon limit log SNR as 7 0. Moreover, it is strange that the capacity 

estimate goes to infinity when 7 tends to zero. Therefore, there are obvious flaws in the result (fTOl i. 

The analytical expression for the conditional probability density function of the channel was obtained in the complex form 
of an infinite series QO), (US, 113 within the Martin-Siggia-Rose formalism based on quantum field theory methods ITSl . In 
the present paper we adopt the Martin-Siggia-Rose formalism and develop a new method for the approximate computation 
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of the conditional probability density function P[y|X]. Using this method we obtain the simple analytical expression for the 
function P[U|X] in the leading and next-to-leading order in the parameter 1/SNR for the intermediate power regime 

N . (11) 

Our method allows us first to derive the analytical expression for the mutual information and then the optimal input signal 
distribution [X] which is different from the half-Gaussian. 

In El a method to calculate the conditional PDF for a nonlinear optical fiber channel with nonzero dispersion in the 
large SNR limit was introduced. Here we illustrate this general approach in application to a simpler nondispersive nonlinear 
optical fiber channel as considered in El, El, El- Since the channel is dispersionless, the temporal signal waveform does 
not change during propagation (note, though, that the signal bandwidth will grow due to the fiber nonlinearity and signal 
modulation). Therefore, instead of considering the evolution of 'ip{z,t) we can consider a set of independent scalar channels 
El, El (per-sample channels) governed by the following model: 

-i'y\'ip{z)\^-ip{z) =r](z), (12) 

where 'ip{z) is the signal function that is assumed to obey the boundary conditions V'(O) = X, ipiL) = Y. The noise r]{z) has 
zero mean {rj{z))n = 0 and a correlation function {'r){z)f]{z'))n = Q5{z—z '), so that the SNR = P/QL, where P and N = QL 
are the per-sample signal power and the per-sample noise power, respectively. The connection between the differential model 
(IT 2 I 1 and the conventional information-theoretic presentation in the form of an explicit input-output probabilistic model and 
appropriate sampling has been discussed in detail in El, m, El- For this per-sample channel we calculate the conditional 
probability density function (in order to illustrate how our method works), the conditional entropy (|5]l, the output signal entropy 
©, and the mutual information (|4]i. Solving a variational problem for the mutual information we find the optimal input signal 
distribution Px[X] maximizing the mutual information in the leading order in 1/SNR. 

The paper is organized as follows. In Section |II] we develop the quasi-classical method for the calculation of the conditional 
PDF P[y|X] for arbitrary nonlinearity in the intermediate power regime (fTTl i in the leading and next-to-leading order in 
1/SNR. We find a simple representation for P[y|X] in this case. This allows us to calculate the output signal distribution 
Pout[Y]. The optimal signal distribution P^*[X] is found in Section |III] Section HVl is focused on the calculation and the 
comparison of the mutual information for various input signal distributions. We demonstrate that there is a range of power P 
where the mutual information IpC^y^x] calculated for a Gaussian distribution P^\x], see Eq. (l34l i below, is closer to Ipopt^^^, 

whereas at large enough power P the mutual information calculated for the half-Gaussian distribution P^\x] is 

closer to /poptr^i than the mutual information We discuss our results in Section [V] 

X I-'I 


II. The conditional PDF P\Y\X] and output signal PDF Pout\Y] at large SNR 
A. ” Quasiclassical” method for the conditional PDF P[U|X] calculation 

The conditional probability density function can be written via the path-integral form El, El, El in a retarded discretiza¬ 
tion scheme, see e.g. Supplemental Materials of Ref. El 


iiiL)=Y 

P[Y\X]= f 

V,(o)=a: 


and can be reduced to the quasi-classical form, see Ref. El: 


■4>{L)=0 


P[Y\X] = e -^ 


Dtp e ^ 


V,(o)=o 


( 13 ) 


( 14 ) 


L 2 

where the effective action 5[-0] = f dz d^tp — i'y\tp\‘^tp , and the function is the ’’classical” solution of the equation 

0 

(5S'['I'ci] = 0, where 6S is the variation of our action The equation (5S'['I'ci] = 0 (Euler-Lagrange equation) has the form 


^ I'fcif = 0, ( 15 ) 

with the boundary conditions = X, = Y. 

In order to find P[U|XJ one should calculate the exponent e Q and the path-integral in Eq. (I141 i. Eirst, we evaluate 
the exponent. To find it we have to calculate the function 4'ci(z) and then the action S'[T'ci(z)]. We found the general solution 
A>ci{z) of ( fTSl l implicitly through the boundary conditions, see Eqs. dbSll-dTSli, and Eq. ( iTTt in Appendix lAl This form of the 
solution is inconvenient for further calculations. Therefore we adopt a different approach and find the solution in the leading 
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and next-to-leading order in 1/SNR, linearizing Eq. ( fTST i in the vicinity of the solution Here 'I'o(^) is the solution of 


the equation (fT^ with zero noise and with the boundary condition '['□(O) = X = \ The function tI>o(z) reads 

T'o(z) = pexp^ip-^ -f , 


(16) 


where p = 'yLp^ = 'yL\X\'^. Note that this solution satisfies only the input boundary condition 'I'o(O) = X = pe'’'^ , and it 

is the solution of Eq. (fTSl l as well. Therefore, we look for the solution of Eq. (fTSl) in the following form 

^cz(-z) = {p + >{{z)^ exp|ip|- -f , (17) 

where the function >c{z) is assumed to be small: jx(z)j -C p- In the general case, the ratio \>c{z)\/p is not necessarily small 
and it depends on the output boundary condition 3<{L). However, the configurations of >c[z) at which 'i'ci{z) significantly 
deviates from '£' 0 ( 2 ) (|>f(z)| ~ p) are statistically irrelevant. Indeed, the expansion S'[4l'o('2^) + i5'I'(z)] oc >c^{z) starts from the 
quadratic term at small xiz), since the action achieves an extremum (the absolute minimum 5'[T'o('Z)] = 0) on the solution 

si'l',,,(»)] 

T'o(z). Thus the exponent e Q and, as a result, the conditional PDE P[F|X] vanishes exponentially if the typical x{z) 
is much greater than ^/QL. 

Substituting Eq. (fTTl i into Eq. (flSl l and retaining only terms linear in x{z)lp, we obtain the following equation which is 
still exact in the non-linearity parameter p: 

(18) 


(Px ^ p dx , p* „ r 1 


The boundary conditions for the function x{z) read 

x{Q) = 0, x{L) = - p = xq + iyo, 

where xq = Re{x{L)} and yo = Im{x{L)}. The solution of the linearized boundary problem (fTSl l. (fT^ reads 

PXQ -yo z ^ (1 - 2p‘^/3)xo + pyo \ z 


(19) 


i?e[xr(z)] = {p 


\ + p^/iL 


Im[x[z)] 


_ ( pxp - yo ( 2p^z 


VI -bpV 3 V 3 L 2 


- 1 


1 -b m73 

(1 - 2pV3)xo 


p- 


■pyo 


1 + /i73 


l) L' 


( 20 ) 


After substitution of the solution Eq. (l20l i in the action we obtain 

{p + x{z)) exp + 7*-^^ I 






dzX — 2iyi?e[>!:] 

1j 


(1 -b 4p2/3)a;2 _ 2pxoyo + i/q 


( 21 ) 


QL(1 + m73) 

Note that here we retain only the terms quadratic in x. However, it is straightforward to calculate the next correction to the 
action (EB which is 0(1/VSNR), see details in Appendix lAl A regular perturbative expansion for x{z) in powers of 1/v/SNR 
can be obtained using the exact equation for the function x{z), see Eq. (l79l l in Appendix lAl 

The next step in evaluation of the conditional probability P[T"|A'] is the calculation of the path-integral in Eq. (fT4b . In order 
to calculate the path-integral in the leading 1/SNR order we retain only quadratic in ip terms in the integrand. Any extra power 
of '(/ or is suppressed by the multiplicative parameter \/QL, because at small Q the main contribution to the path-integral 
comes from ip ~ \/QL. Moreover, as soon as we calculate the path-integral in the leading order in Q, we can substitute 'I'o(z) 
for 'I'c/(-2:) in the action difference 5'['I'ci(2:) + '0(-z)] — S'[4'cz(z)]. To find P[y|Ai] in the next-to-leading order in 1/SNR we 
should retain both x{z) in 'i'ci{z) and higher powers of ip in the action difference. Details of the path-integral calculation in 
the leading and next-to-leading order in 1/VSNR are presented in Appendix |B] Taking into account the expression for the 
action dSSl l and the result of the path-integral calculation (I1091 l we obtain the final result for P[y|2f]: 

f {l + 4:p‘^/3)x'^-2pxoyo + y^ \ 

PiY\X] = -i + m73) -£ I^ (p(15 + p^)xo - 2(5 - pV^)yo) - 


ttQLa/I + M^/3 


15(l-bpV3)2 


p/p 


135QL(1 -b m73)' 


3 ( 5 /-b 33/4^-b 30) 


-|p (V + ISft^ + 225) -b ( 23 / -b 255p^ - 90) + P (20/ -b 117p^ - 45) xoyo 


O 


( QL 


O (773q|X|2) 


( 22 ) 
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where xq and yo are the functions of X and Y defined in (fT9T l. Since we consider here only the result (l22l i for the large 
SNR limit we imply that |Xp ^ QL. Note that the conditional PDF P[F|X] was already derived in ifTSi in the form of an 
infinite series. Our result (l22T i for the function P[F|X] is the analytic summation of this series in the limit of large SNR and 
intermediate power region 

QL<P«C (QLV)”^- (23) 

The accuracy O in Eq. (l22l i appears in the calculation of the path-integral as a result of neglecting higher powers 

of the field ip, see Eq. (|90] | in Appendix iBl One can show that the normalization condition J DYP[Y\X] = 1 is fulfilled. Also 
one can check that the distribution (|2^ obeys the following important property 

nm^p[y|x] = (24) 


The expression (l24l i is nothing else, but the deterministic limit of P[y|X] in the absence of noise. Also Eq. (|2^ has the 
correct limit for the linear channel (7 —0); 


p(°)[y|A:] 


^-\Y-Xf/QL 

ttQL 


(25) 


that is nothing else but the conditional PDE for the linear nondispersive channel with AWGN. 


B. Output signal PDF Pout\X\ 

Now we proceed to calculate the probability density function of the output signal Pout\Y]. Let us consider the integral, see 
Eq. 0, 

Pout[Y] = J VXP[Y\X]Px[X], (26) 

where the function Px [2f] is a smooth function with a scale of variation P which is much greater than QL. In that case we 
can calculate the integral ( l26b with accuracies O (1/SNR) and O (j^L^QP^ by Laplace’s method ll20ll . see Appendix ICl The 
result has the form: 


Pout[Y] = / VXP[Y\X]Px[X] = Px 




(1 + 0 (1/SNR) + O {-f^L^QP)) . 


(27) 


The main term of this result (l27l i can be easily obtained from the following reasoning. The function P[y|X], see Eq. (|2^ . 
varies on a scale of order QL which is much less than the scale of Px\X\ (the function P[y|Ai] is essentially narrower than 
the function Px[X]). Also P[y|A'] has the delta-function limit (l24li and therefore in the integral (|2^ it can be replaced with 
the delta-function. Note that to obtain the result (l27l) we do not require the limit Q —0 but only the relation between the 
scales P and QL to be satisfied. In what follows we will omit the terms 0{. ..) for brevity and restore them in the final results. 
Eor the case of the distribution Px[X] which depends only on |Ar| we have Pout\Y] = Pj(:[|y|]. Lor such distributions we can 
calculate corrections to ( l27b in the parameter QL in any order in QL. 

Let us restrict our consideration in the remainder of this Section to the distributions Px[X] depending only on |Ai|. We can 
use the P[y|Ar] found in Ref. ifTSll . see Eqs. (11)-(13) therein. In this case Pout\Y] is a function of |y| = p' 

fdppe ,28, 

0 


where Io{z) is the modified Bessel function of the first kind. Using this representation we can obtain the simple relation for 
Pout[p'] calculation in the perturbation theory in QL. To this end we perform the zero order Hankel transformation ll20l : 


P[k]= J dppJo{kp)Px[p]. (29) 

0 

of both sides of Eq. ( |28] |. then we use the standard integral ED with Bessel and modified Bessel functions 


dzze Jv (bz)P(cz) = —Ju 

2p 


he \ 


-bd 

I e '‘p 


and arrive at the simple relation between the Hankel images 

Pout[k]=e-'^"^P[k]. 


( 30 ) 
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Performing the inverse Hankel transformation 


Px[p] = J dkkJo{kp)P[k], 
0 


we obtain 


( 31 ) 


(32) 


where Ap = 


Pout[p] = e 4 ^pPxlp], 

_ __ ____ _ _ _radial Laplace operator. From the relation ( l32b the problem of finding (QL)" 

corrections to Poutfp] reduces to the exponent expansion and straightforward calculations of the action of the differential 
operator A" on Px[p]- 

Let us consider the widely used example of the modified Gaussian distribution 

exp{-Pp'^/{2P)} pf^-'^ 


Pflp] = 


(33) 


ttT (/3/2) (2P//3)'^/^ 

For f3 > 0 the distribution [p] is normalized to unity, 27r dppP^^ [p\ = 1, and has the average power P, 2tt dpp^P^^ [p] 
P. The distribution P^\x] generalizes the half-Gaussian distribution ® for /? = 1 and the Gaussian for fi = 2: 




-|xF/p 


(34) 


Inserting 0^ into Eq. 
has the form: 


we obtain a standard integral which can be found in 11211 . The result for the output signal PDF 


.(/3)r 


/3 


\Y\^2P 


out[Y] (^2’^’QL(2P + ^QL) 

exp{-|y|VQi} ( PQL 

2P + I3Ql) ’ 


ttQL 


(35) 


where iF'i(A 1] z) is the confluent hypergeometric function that reduces to for the Gaussian case and to e^/'^Io{zj2) for 


the half-Gaussian case: 


p!>11{y] = 


7r{P + QL) 


exp< - 


|r|^ 


P + QL 


(36) 


_ exp{-|y|V(2P + QL)} ^ 
7t\Y\,/7t{2P + QL) 


'Trlyp 

_!-!_g QH2P + QL) /g 

QL 


\Y\^P 


QL{2P + QL) 


(37) 


Note that the result for in Ref. ifTH . see Eq. (38) therein, for the half-Gaussian distribution is incorrect. To demonstrate 

the general result of Eq. (l27l let us consider Eq. (l37l in the case QL <C \Y\'^ ~ P. Eor the case one can obtain: 

/ 1 


The result (l38l) coincides with Eq. (l27l) with the accuracy O (j'^L^QP). 


3 ( 1 ) r 


(38) 


III. Optimal input signal distribution at large SNR 

The optimal input signal distribution at large SNR can be found calculating the mutual information (|4|i and then maximizing 
the result with respect to the input signal distribution function Px [A]. Let us start from the calculation of the output signal 
entropy L[[Y], see Eq. (|6]l, at large SNR. When the parameter SNR ^ 1 we can substitute Px [Y exp {—i 7 |FpL}] instead 
of Pout[Y] due to the relation ( |27] |: 

277 00 

H[Y] = - j d(t>j dp'p'Px [p'P*] X 


0 0 
\ogPx [p'P*] +0 


1 


SNR 


O {pL^QP) . 


(39) 
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In order to obtain Eq. ( l39l l we have performed the change of the integration variable <j) = + 'y\Y\'^L. One can see that 

the output signal entropy ( |39] ) coincides with the input signal entropy in the leading order in 1/SNR and j^L^QP. 

The conditional entropy iJ[y|X] can be calculated by substitution of P[y|X] in the form of Eq. (l22l i into Eq. ©. After 
the substitution we change the integration variables VY = ctReYdlmY to dxodyo. Then we perform integration over xq, yo 
and obtain 


H[Y\X] = 1 + log(7rQL) + 


271- 




dppPx 


pe 


i4>' 


(X) 


log ( 1 


rL^ 4 
-^p 


o 


SNR 


O {-f^L^QP) 


(40) 


0 0 

where the first two terms come from the Gaussian type integrals over VY in the conditional entropy definition Q and the 
normalization factor ttQL in Eq. (l22l i. The third term in Eq. ( l40l i comes from the normalization factor \/l + At^/3, see Eq. 
(I 22 I 1 . Note that there are no terms which are O ( l/v/SNR J in Eqs. ( |39] | and (l40l i. Indeed, the integrals with the odd powers 
of xo and yo vanish when integrating over xq, yo in Eq. @ for 

To find the optimal distribution P^*[X] normalized to unity and having a fixed average power P one should solve the 
variational problem for the functional J[Px, Ai, A 2 ] 


J[PxAiA2]=H[Y]-H[Y\X]-Xi (^J VXPx[X]-1^ - 

X 2 (^J VXPx[X]\X\'^-Py (41) 

where Ai ,2 are Lagrange multipliers. We substitute H[Y] and iJ[y|X] from Eqs. ( |39] | and (l40l) to dril l, perform the variation 
of the functional J[Px, Xi, X 2 ] over Px[X], Ai, A 2 , and write the Euler-Lagrange equations 6 J[Px,Xi,X 2 ] = 0: 


j VXPx[X] = 1, 

J PXPx[X]|X|2 = P, 

-1 - logPx[X] - i log (^1 + - Ai - A2|X|2 = 0 . 


The solution P^*[X] of Eqs. (I42l) - (l44l) referred to as the “optimal” distribution depends only on \X\ and has the form: 


PTiX] = iVo(P) 


exp{-Ao(P)|X|2} 


where functions Ao(P) and Ao(P) are determined from the conditions (l42l i. (|4^ : 

o 

dp 

\/1 + 72LV/3 


J PXP7‘[X] = 2-kNo{P) j 


= 1 , 


(42) 

(43) 

(44) 

(45) 

(46) 


yPXPj‘[X]|Xp= 2ttNo{P) j- 
In a parametric form this dependance reads 


dp pY 


3p-MP)P 


7P 


y/r+T^IV/S 

7P 


= p. 


nV^Giay 


(47) 


(48) 


here G{a) = = f {Po (a) — To (a) I with To(q:) Ho{a) being the Neumann and Struve functions of zero 

order, respectively. The parameter a(P) > 0 emerges as the real solution of the nonlinear equation ^ log G{a) = —7PP/v/3, 
which comes from Eqs. (l46l l and (l47l i. Let us emphasize that the optimal distribution obtained here, P^*[X] (|45]) . is different 
from the half-Gaussian distribution, see Eq. (l33]) for /3 = 1, whereas in the Ref. IH the half-Gaussian distribution was 
considered as optimal. Eor sufficiently large values of the power P, such that log{'jPL) ^ 1, we can simplify (l48l l using the 
asymptotic expansions of Yq (a) and Pq (ct) at small a, see ED: 


Ao(P) 

No{P) 


1 - loglog(C' 7 )/log(C' 7 ) 


7^1og(G7) 


-Ao(P), 


(49) 



















where C = 2e and 7 = 'yLP/ v^. At small P, the parameter 7 ^ 1 , the solution of the Eqs. (l46T l and (l47l i has the form: 

HP) = p (1 - 27 ") , No{P) = ^ (1 - f) . (50) 

It is worth noting that at 7 —0 our distribution (|45]) approaches the Gaussian distribution (l74l) that is known to be optimal 
for the linear channel IT]. 


IV. The mutual information 

In this Section we present the entropies and the mutual information for P^\x\ and P^^[X] and compare our new results 
with those already known. 

Substituting the expression (l45t for P'^^\X] in equations (I39b- (l40l) and using the definition (HI we obtain the mutual 
information for the optimal distribution in the leading order in 1/SNR: 

Iprm = PHP) - - log(7reQL) + O + O {i^L^P) , (51) 

which gives a capacity estimate in a wide range (l23l l of the average power P. The mutual information (ISTT i is depicted 
by the black solid line in Fig. [T] as a function of power P for the following parameters: Q = 1.5 x 10“^mWkm“^, 
7 = 10“^ mW~^km“^, L = 1000 km. For these realistic parameters, the power range (l23l l is actually very wide: 

1.5 X 10"^mW < P < 0.66 x lO^mW. (52) 




Figure 1. The mutual information for various input PDFs as a function of input average power P for the parameters Q = 1.5 x 10 ' mWkm 
7 = 10“^ mW~^km“^, L = 1000km. (a): The solid black line, blue dashed line, red dashed dotted line con'espond to the optimal PDF P^^[X], 


Gaussian PDF P^\x], and half-Gaussian PDF Py^[X], respectively, (b): The solid black line con'esponds to I 


Mx]’ 


: Eq. (5T); the red dashed dotted 


line corresponds to the mutual information for the half-Gaussian distribution / (i) , see Eq. {57) for /3 = 1; the red dashed horizontal line coiTesponds to 

our limit (ID at 7 ^ 1 for the half-Gaussian distribution; the black dotted horizontal line corresponds to the result □3, see Eq. m- 


There is no simple analytical form for Nq{P) and Xq{P) therefore to plot Fig. [T] and Fig |2] (see below) we calculated 
Ao(P) and Nq{P) numerically. For large and small values of the parameter 7 we can use the solutions in Fqs. (Hill and (|50ll, 
respectively. At small 7 = jLP/\/3 we obtain 

Iprix] = log (1 + SNR) - f + OH) + O , (53) 

which is simply the Shannon capacity, log (1 + SNR), of the linear AWGN channel ([T]) with the first nonlinear correction. 
In Eq. (I 53 I 1 the unity in the logarithm is beyond the accuracy of our calculation but we keep it to bring to notice that the 
derived expressions dTOl i and (l53l l have the correct limit when the parameter 7 tends to zero (in contrast to the Eq. (35) in 
Ref. ifTTIl l. In Eq. (fsTI) we omit the accuracy O {pi^LP'QP) since the parameter 'j^L^QP is of order of 7 ^/SNR. In the second 
power sub-interval ( 7 ^)“^ <C P ’ using Eq. ( 1491 ) one can see that the mutual information increases very slowly 

(loglog) with P 


‘[X] = - log {QL^i) - 1 + 
1 


logs 


loglog 


C'7PP\ 

Vs ) 


log (CjLP/H) 
0{l/\og^{jLP)) 


loglog ( £27) + 1 _ loEbgiCLP/^) 

^ Vs y log C7PP/V^ 


0{l/SNR) + 0{j^L^QP) , 


(54) 



















9 


as opposed to the constant behavior of the mutual information for Gaussian-like distributions of an input signal, see formulae 
( |58] | and ( |59] | below. 

In the remainder of this Section we perform an analysis of the mutual information for the distribution P^\x], see Eq. (l3^ . 
generalizing the half-Gaussian distribution (O (see, for example Ref. lfT4ll ') and the Gaussian input PDF 0^ . In the leading 
order in 1 /SNR from ( l39l l we obtain 


H^[Y]= log 





2-P 

2 



(55) 


where '(/'(z) is the digamma function = r'( 2 :) /T{z), where ■0(1) = —1 e and -0(1/2) = —^e — 2 log(2). The substitution 
of Eq. (|33 T i into Eq. (l40l i gives 


Hp[Y\X] = log (ttcQe) + 

oo 


(56) 


with at least 0(1/SNR) accuracy. The integral in Eq. ( |56] | can be calculated analytically using Ref. ||2T1 . however, the result 
of the integration is cumbersome, hence we do not present it here. One can easily obtain the mutual information by 

subtracting Eq. (l56T l from Eq. (I55b : 


/pW) [X] = log SNR + log _ 

OO 

^(l-^(f))+o(^)+O(.V0P). (57) 

The mutual information is depicted in Fig. [Ha) for the Gaussian distribution by the blue dashed line, and for the half-Gaussian 
by the red dashed dotted line. One can see that at small P the mutual information for the Gaussian distribution is greater than 
that of the half-Gaussian, whereas at P > llmW the mutual information is greater for the half-Gaussian distribution. Note 
that is greater than for all values of P, as it should be. At 7 ^ 1 the mutual information IpW^^x] t^^es the 

form 


r 1 \ 2 - /3 logs , (I3\ , 

log(r(/3/2))+o(^^)+0(72 l3QP). (58) 

One can see that at large SNR goes to a constant in the interval of power P considered, and this constant depends on 

the noise power QL. We remind tiiat Ipopt increases as log log P in the region under consideration. The mutual information 
for the half-Gaussian distribution ® in the regime -7 ^ 1 can be obtained as a particular case of (l58l l for /3 = 1; 


[X] = “ {QL^l) + log 2 -f 

^ (^) + ^ ■ ( 59 ) 

Comparing our expression (|59] | with the result (40) of Ref. m we have an extra term + log 2 due to our more accurate 

calculation of iJ[y|X]. Our result (l59l l and the result of Ref. Ifl4l . see Eq. (fTOb . are presented in Fig. [Ub). In Fig. [Ttb) one 

can see that the mutual information dSTT i for the optimal distribution exceeds the limit ( |59] | at P ~ 190mW. At this power the 

difference between the limit ( l59b and evaluated on the base of Eq. (l57l) with /3 = 1 is of order of 1.5% and getting 

smaller at higher P. 

Since we have now found P^*[X] in the power region (|2^ . we can calculate an approximation for the capacity of the 
considered per-sample nonlinear channel. By definition it coincides with the mutual information expression (ISTl i: 


C — Ipopt . (60) 

Let us emphasize that this result for the capacity has accuracy O (1/SNR) -f O {^'^Li^QP'), see Eq. ( ISTT i. The comparison of 
the approximation (|60] | with the Shannon capacity of the linear AWGN channel is presented in Fig. |2] One can see that the 
Shannon capacity of the linear AWGN channel is always greater than the approximation dhOl l for the nondispersive nonlinear 
fiber channel for the considered region of P. 
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Fiaiirp 9 SIhannnn ranaritv anH thp ranaritv nf thp nnniinpar rhannpl T for the parameters Q = 1.5 X 10 ^ mW km 7 = 10 ^ mW ^km ^ 

, the red dashed-dotted line corresponds to the Shannon limit log[l -|- SNR], 


V. Conclusion 

We have developed a new approach to the calculation of the conditional PDF via the path-integral representation (fT4l) at large 
SNR for the intermediate power region (|2^ . This may be an especially useful technique for complex nonlinear channels in 
which the calculation of the conditional PDF is technically challenging. Applying our method to the per-sample nondispersive 
nonlinear fiber channel, we derived compact analytical expressions for the conditional PDF, conditional entropy and the 
entropy of the output signal for different input signal PDFs Px[X]. Moreover, we solved the variational problem on Px[X] 
maximizing the mutual information in the leading order in 1/SNR in the power region (|2^ . That allows us to find the optimal 
input signal distribution (|45]) and the approximation for the channel capacity (BTI) up to 1/SNR and ^‘^L^QP corrections in 
the power interval QL ^ P <C which is extremely wide for realistic parameters, see (l52]) . The found distribution 

P‘^*'[X] is different from the half-Gaussian one, and at the zero nonlinearity P’^*'\X] approaches the Gaussian distribution. 
We demonstrated that the approximation found for the capacity of the channel considered here dbOl) is always greater than the 
mutual information calculated for Gaussian and half-Gaussian distributions, and lower than the Shannon capacity of the linear 
AWGN channel. 
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Appendix 

A. The classical solution T'cz and the action S'[T'ci]. 

In Ref. ifTTll we have shown that in the case SNR = P/QL 1 the conditional probability can be written in the form: 

^{ L)=0 


P[F|-A] = exp < — 


Q 


J Vij} exp \ — 
^(o)=o 


S'[T'ci( 2 ;) -I- i>{z)] - S'[T'c/(z)] 

Q 


where for the nondispersive model the effective action reads 

L 


= j dz 


dzTp - 


The action (162b is associated with the l.h.s. of the nonlinear Shrodinger equation 

dzip{z) - i-f\ip{z)\‘^ilj{z) = r]{z), 

where the noise r]{z) has the Gaussian nature: 

{r]{z))r, = 0 , {r]{z)fj{z'))r, = QS(z - z'). 

The measure in Eq. dMb is defined as 

. N-l 


'D'ih = lim (— —— TT dReipidlmibi^ 
A->o V AirQ/ -*■ -*r 


(61) 


(62) 

(63) 

(64) 


( 65 ) 
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here ipi = tp{zi) and A = ^ is the grid space. 

Now we consider the difference of actions in the exponent of the path-integral in Eq. dMl . 

+ ' 0(^)1 - = 

L 

y 10^0 - i 7 ( 20 |^'ciP + + 27/TO [{dz'^ci + * 7 ^cz|^'dP) (2^'ci|0|^ -f + 

0 

7^ |2^'c;|0p + +'0|'0n^ -f 27 /TO i7(20|«'ciP + 0’j'ci)) (2^'ci|0p + +0101^)] 

27 /TO [[d^^ci + *70cj|^c/p) 0 I 0 P] |- 


( 66 ) 


In Eq. dMT) the function '^ci{z) is the solution of the equation (5iS'['I'cj] = 0 (Euler-Lagrange equation) which has the form 

= (67) 




dz^ ' I ‘-‘I 

with boundary conditions 'l'ci(O) = X = |A| exp[i0(^)], 4'cz(i) = Y = |F| exp[z0(^)]. It is easy to find the solution of 
Eq. (l67l i in the polar coordinate system: 'l'ci(-z) = C = z/L. The solution depends on four real integration constants. 

We denote them as E, jl, (^0 and Oq. There are two different regimes of the solution: in the trigonometric regime one has 
E = ^ > 0, and in the hyperbolic regime E = —^ <0. Eor both cases instead of E we introduce the non-negative parameter 

fc= V^- 

,2 

• In the trigonometric case (E = ^ > 0) we have the solution for p > fc > 0: 


P^(C) = ^ (a + cos[2k{C - 


0(C) = f (C - Co) + 


sin[2fc(C- Co)] 
4/ 


arctan 




Here the integration constants /i, k and Co must be found from the boundary conditions: 

|Xp = p2(0) = (^j 2 + y/fi^ -k^ cos[2A:Co]) , 

^ ■ (^0 -f v//i2 - cos[2/c(l - Co)]) , 


\Y\^=p^il) = 


2 L 7 ' 

d-, 


c,(^) = 0(0) = -§Co - 


sin[2fcCo 


— arctan 


C,(^) = 0(1) = |(i - Co) + v/0^^ 

Then one can find the action 


Ak 

sin[2fc(l - Co)] 
4fc 


( 0 - - A: 2 ) 


tan[fcCo 


arctan 


(0- \/02 - fc2) 


+ 00 j 

tan[A:(l - Co)] 




S[A!ci{z]E= Y,/i,Co,0o)] = ^ ( 0 - \/p-k‘^ 


sin[2fe(l - Co)] + sin[2fcCo] 

2fc 


( 68 ) 

(69) 

(70) 

(71) 
+ 00- (72) 

(73) 


In the hyperbolic case (E = — 0) we have the solution for fc > 0 and arbitrary jl in the following form 


7 (c) = ^( - + 77Tr^cosh[2C(i: - Co)]), 


C(C) = -f(C-Co) + 7A^^^®fA!)l 


— arctan 


{ji + ^/WT^) 


tanh[fc(C - Co)] 


(74) 


00 , 


where jl, k, Co, and Oq are derived from the same procedure as in the trigonometric regime. The action reads 


fc 2 




S[Aiciiz;E = —^,A,Co,0o)] = 2 y 




I p 2 sinh[2fc(l - Co)] + sinh[2fcCo] '\ 


(75) 


Note, there are two solutions of Eq. dFTl) with constant p{z) = p(0) = p obeying only the input boundary condition 4'o(0) = X. 
The first one reads 

(76) 


T'o( 2:) = pexp + * 


where p = 'yLp^ = ^L\X\'^. This function corresponds to the solution representation dhSl) with fc = 0 and jl = p or Xo the 
solution representation (l74l) with k = 0 and jl = —p. The function 4'o(2]) is the solution of the Eq. (l 6 ^ with zero noise and 
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with the input boundary condition. Furthermore, ^'o(z) delivers the absolute minimum of the action (l62l i: S'[4'o(z)] = 0. The 
second solution of Eq. ( |67] | with constant p{z) is the trigonometric regime dhSl ) case with jl = k = 2p: 

'^p=const{z) = , p = ^Lp^ = jL\X\'^. (77) 

To find the solution of Eq. WT\ one should express the integration constant through the boundary conditions. Instead, we 
exploit the fact that the noise power QL is much less than the input signal power P = = p^. In other words, we will 

find a solution of Eq. d&Tl i that is close to 'I'o(2): it is the solution of Eq. (l6^ with zero noise which provides the absolute 
minimum of the action 5[T'o{2)] = 0. In that fashion we perform the substitution in Eq. d&Tl i: 

^ci{.z) = (p + >ir( 2 ;)) exp , (78) 

where the function x{z) is assumed to be small: >c{z) p for all configurations of 'i!ci{z) providing S[^ci{z)\/Q = Oil) 
when QL tends to zero. We have the following equation on >c{z) resulting from the Eq. (l67] i: 


dpK a d>c 
dz^ dz 


. ft „ . , 


+ lOixp + 3>r2] + 


L2p2 


, T dx 

AiL ——h 9px + lApx 
dz 




■x-^ + 


L2p3l 


[3|; 


2xr^ 


3/j,^ 

Zv 


\x\'^x. 


(79) 


We present x{z) as a perturbation theory decomposition in powers of l/VSNR: x{z) = xi{z) + X 2 {z) + ..., where xi{z) 
is of l/v/SNR order and provides the leading order contribution, X 2 {z) is of 1/SNR order, and so on. 

• The linearized equation for the function xi(z) = xi(z) + iyi(z) can be obtained from Eq. ( l79b by omitting the r.h.s. of 
this equation: 


d^xi „ u dxi , „ r 1 


(80) 


The boundary conditions T'cz(O) = X and 'l'ci(l) = V = lead to 




p- 


, dsTl) is polynomial 


>^i(0) =0, 

xi(L) = xo + iyo = p'e*' 

The solution xi{z) = xi{z) + iyi{z) of the linearized boundary problem 

2 : 1 ( 2 :) = ( - pai{X,Y)^ +a 2 iX,Y)^^, 

l(2:) = ( — + pa2iX,Y)— + ai(X,F)^—, 


yi[ 


where coefficients ai{X,Y) and a 2 {X,Y) can be found from the boundary conditions (ISTT l and have the form: 


ai{X,Y) = 


-pxo + yo 


a 2 (X,Y) = 


(1 - 2p‘^/3)xo + pyo 


l + k/3’ l + p2/3 

with xq = Xo(X^ Y) and yo = yo{X^ Y) being determined from Eq. dSTI ). In the leading order the action reads 

L 


Q 


px) 


A>o{z) + xiiz)P>^i+^‘^ 

(1 + 4p^/3)af — 2paia2 + a| 

gZ 




dzXi — 2iYRe[xi] 

1j 


^4 


+ 0 


f 1 


VksM 


(1 +4^73)2:0-2^2:02/0 +2/0 , q ( 1 


QL(l+/i2/3) 


VVsM 


(81) 


(82) 


(83) 


(84) 


Let us proceed to the next-to-leading order correction to P[F|W]. We should calculate the next approximation >!: 2 (- 2 :) to 
the solution (l78l l. Taking into account Eq. dSOl l we present the equation for X 2 iz) in the form 

d?X 2 


. p dx 2 


14 ^ ~ ■ 4f2^ek2] - 42^ + 3 ^,) ^ ^ [5k + lOkiP + 3^1^] 


(85) 
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where the boundary conditions for >C 2 {z) read >r 2 ( 0 ) = X 2 {L) = 0. The solution 3 <2{z) = X 2 {z) + iy 2 {z) of Eq. ( l85T l is 
polynomial in z and quadratic in xq and uq: 

270(1 + /r2/3)3l^ l ) 

|/r (2/ - 15p2 + 585) + 2 (iS/i^ 15 ^ _ igq) xoyo + y. ( 2 ^^ + 15 ) ( 5^2 _ 9 ^ y 2 _ 

5 + 3) I; (m - 15) aig - 4 [y^ - 6) xoyo + M + 9) j/q) + 

.2 (86) 

by {y^ + 3 ) ^ (3 (5p^ - 3 ) aig - 36yxoyo - - 15) y^) + 

z3 

20/r^ [y^ + 3) j^{yo - y-xo) {“^yya - {y^ - 3) xq) - 

z^ ^ 

20 y^ {y^ + 3 ) jj{yo - yxof^ 


( \ p/p 2:N Z 

270(1 + mV3)3 V l) 

{ (7/i^ - 75/i2 + 360) xl + 6/i (/r^ + 75) x^y^ + 3/i2 {by'^ + 39) t/g+ 

2y {{y^ - 4 / + 255p2 + 180) xl +y {y^ + 15) {liy^ + 3) xoyo + (5/ + 36/1^ - 9) yl) - 

1j 

z2 

14p {y^ + 3 ) j^iVo - yxo) ((15 - Ay^) Xq + 9yyo) + 

3 

( m ^ + 3 ) ^{yo - yxo^y 


In the leading, see Eq. ( l84l i. and next-to-leading order in l/VSNR the action reads 


Q 




(1 -f 4pV3)xg - 2yxoyo + yl 
QL(l + p2/3) 


P/P 


135QL(1 -f/r2/3)' 


-{Ai(4/-f 15/i2-f 225) xg-f 


(23/ -f 2bby^ - 90) /yo + P (20/ -f 117/ - 45) xoyl - 3 (5/ -f 33/ -f 30) /| -f O 


( 1 


VSNR 


(87) 


( 88 ) 


B. The path-integral calculation. 

To calculate the conditional probability density P[F|X] in Eq. dMl l one should find the pre-exponent path-integral, referred 
to as the quantum corrections near the classical solution '^ci{z), in the leading and next-to-leading order in 1/VSNR: 


i/.(i)=0 

Iqc['^ci{z)] = J Vtj} exp {- 
y’(o)=o 


5[T'e;(z) -f '0(z)] - 5'[T'ei(z)] 

Q 


(89) 


In what follows we are interested in the leading and next-to-leading order corrections for the path-integral (IhTT i. That is why we 
retain only quadratic in ip terms in Eq. ( l66l l. All these terms are placed in the second line of Eq. ( l66l l. As it will be demonstrated 
below an extra power of ip results in an extra power of ^QL. In the leading and next-to-leading order calculation of the path- 
integral we should take into account the first correction (>ri(z) oc y/QL) to the solution 'l'ci(z), see Eqs. (fTSI l and (l82] |. Now we 
put (l78l l with >ri(z) and ip{z) in the form ip{z) = u{z) exp {iyj^ + into the first line of Eq. (|66]l. In our approximation 

we obtain 


L, 

S'[T'c;( 2;)-f'0(-2)] - <S'[T'c/(-2)] = y^2;| - i^(u-f u) 


2—Jm 


Lp 


2 (j)zU -j- i^iu u 


■Pi 


t(>ri -I- >ri) -I- {^dzXi + *^(^1 + 


2|ui 


O 


+ O {Ipu^>c2) + O (//n^) }. 


(90) 














14 


We substitute this difference in the exponent in Eq. (l89l l. Then we expand the exponent at small Q and obtain: 

L 


exp< - 




Q 


^ = exp ^ —— J dz dzU — i^{u + u) || 


1 - 


QLp 

O 


Tv 

J dz 2^1 
0 

^)+Oh-0LV)}. 


dzU + + ^)) (+ ^i) + 


) + (a 


>ri + *^(^1 + ^i)) (2|u|^ + 


(91) 


Here we imply that any extra power of u or x: is suppressed by the multiplicative parameter \fQL, because at small Q the 
main contribution to the path-integral comes from u ^ \/QL. We substitute this expansion (IMT l into the path-integral ( |89] | and 
change the variable from '4>{z) to u{z) and arrive at 

ii(L)—0 L 

OzU — l—{u + u) J 


Iqc[^o{z)] = y exp I - — 


dz 


1 - 


!i(0)=0 


2fi 

QLp 


Tv 

Im J dz 2(j)zU + (^u{>ci + xri) + uki^ + (j)z^i + 


2|uP -I- 




(92) 

To calculate the leading and next-to-leading order contributions to /qc['I'o(^)] in 1/SNR we should take the first and the 
second terms in the square brackets in Eq. (|92] |. respectively. We start our consideration from the leading order. In this case 

in the retarded discretization 


we represent the path-integral /q(/['I'o(-z)] = Ju(o)=o ~ Q Jo 

scheme: 

u{L)=0 L 


AO) 

‘qc 


{-^lo 

[T'o(2)] = J Vu exp^- ^ J dz dzU- i^{u + u) | 


dzU - 


(u-|-u)| I i 


lim 


N —^oo \^7tC^L 


N 


m (0)=0 
N tV-l 


r ( AT ^-1 ) 

I n {“I: S 1 ' 


where we use the measure (l65] t and the notations u{zj) = Zi = Ai, A = ^ and = 0. 


'^3 ' "^^3 

(2) (2) (2) 

The sequential integration over 'w)v- 2 ’ • ■ Wi is trivial: 




[A-Bf 
2(ti -I- T 2 ) J 


It leads to the remaining integral (over uA\ i = 1,N — 1) of the form 


lim 
JV — >^00 \^ 7 ^(^L 


N (ttQL/N)^ 


N~1 ^ N — 1 




N 


N-1 


Vn 


W: 


,(i) 


Z=1 


QL 


(94) 


(95) 


jj'=i 


where we denote a = , and the (iV — 1) by {N — 1) matrix M{a) has the following elements: Mi^i = 2 -f a, 

Mi^i±i = —1 + a, i = 1,..., N — 1, Mij = a, j ^ i, j ^ i ± 1. It is straightforward to calculate the determinant of M{a) 


and hence to perform the Gaussian integration over u\ 


( 1 ) 


det[M(a)] = N + a 


N^{N^ - 1) 


12 


(96) 


4c['I^o(^)] = 


■nQL^Jl + pP'/2> 


(97) 
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To calculate the next-to-leading order contribution to the path-integral ( |92] ) we should take the second term in the square 
brackets in Eq. ( |92] |. To find this correction we should calculate the integral (the correlator): 


{u^^\z)u^^\z')) = 




iSkMz)], 


s 

II 

0 



dzU—i^(u-\-u) 


u^°‘\z)u^f^\z') = QLG°‘'^{z,z'), 


(98) 


(o)=o 


where we have introduced the dimensionless Green matrix G°‘'^{z, z'), a, /? = 1, 2. The standard method for the Green matrix 
calculation is the calculation of the generating functional IHl: 

u(L)=0 L L 

Z[Ji,J 2 ]= Jvuexp^-^J dz dzU - i^{u + u) + J dz (^Ji{z)u^^\z) + J 2 {z)u^‘^\z)^, (99) 

u ( 0)=0 0 0 

then any correlator can be derived from the variation of the ^[Ji, J 2 ] over Ja, for example 


{u<^‘^\z)u^P>{z')) = 


WfZW = 


1 


SZ[Ji, J2] 


J 2 —0 


( 100 ) 


^[Ji,J 2 ] 5J^{z)5Jp{z') 

The calculation of the generating functional can be performed in the same way as the calculation of the normalization integral 
(|9^ : the integration over u'^'^ followed by the integration over ■ The only new element in the calculation of the Gaussian 
integrals with the sources is the inverse matrix M{a)~j for M{a)ij = a + 2Sij — Sij+i — Si+ij defined herein above, 
see the text after Eq. ( |95] |. The calculation is simple (after the observation that det[M[a)]M{a)~^ is linear in a), and we 
only set out the result 


Mia)-} = N 


N n) ^ N 4det[M(a)]fV 


* Afi_A 

N N \ N 


,( 101 ) 


where det[M(a)] is given by Eq. ( |96] |. and limAr_>.oo (Q:iV"‘/(4det[M(a)])) = 3/r^/(3 -I- We present the result of the 
generating functional calculation in the form of a Green matrix convolution with the sources Jc,: 


Z[Jl,J2] = 


1 




exp 



dz j dz Ja{z)G°‘’^{z, z)Ji 3 iz) 


( 102 ) 


0 0 

where the Green matrix is Hermitian and it has the following elements: 


G^^\z,z') = G^^Hz\z) = Biz' - z)_ 1 - - - 


2L 


3p^ 


LJ 4(3+ ^2) 




G^'‘^iz,z') = G‘^'^iz',z) = 


2(3-b/i2) 


1-1 


L L \L ) 


Hi-') 


z 

IL 


-b 


(104) 


z') = = 


Z\z' ^ „ 9 /, ^zz' 




Biz-z') 

6(3 4- 

4- {z O z'}. 


(105) 


The second way to obtain the expression for the correlator ( l98l l and Eqs. (I103b - (ll05b reflects the fact that the Gaussian 
integral ( l99b is saturated in the vicinity of the saddle-point solution of the equation of motion (i.e. the Euler-Lagrange equation 
for the action in question) im. Thus to find it we should solve the set of equations 


K^^Gz^\z,z') = jbibiz-z'\ 

where the matrix differentiation operator K for the functions uiz = 0 ) = u(z = L) = 0 is defined as 


L, 

dz dzu - i^{u + u) = J dzu^°‘'^iz)ka, 0 U^'^\z), 


(106) 


(107) 
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and it has the form 


K = 2 




-a? 


(108) 


The boundary conditions for equations (IIO 6 I 1 are as follows; = 0 , 2 :') = G°‘’^{z = L,z') = 0. The problem has 

the unambiguous solution (I1031 l- (ll05b . Note that the homogeneous solution of the Eq. (IIO 6 I 1 is governed by the solutions of 
Eq. dSOb obtained above. 

Using the correlator dMI l with (I103b - (ll051 l one can easily calculate the first correction presented in the second line of 
Eq. ( l92b . This term is proportional to >ci{z) oc y^QL hence delivering the first correction in l/v/SNR to the leading term 
(l97b . The subsequent integration of the elements (I103b - (ll05b with the solution (l82t for }<i{z) is trivial, however the proper 
way to understand the discontinuous derivatives of the Green matrix elements (I103b - (ll05b at the same point z' = z is the 

dzG°‘’^{z + 0, z')\z’=z- Finally we have 


retarded scheme adopted in our approach E): dzG^’P{z,z')\z'^z 




O 


QL 

n2 


QLv'I + mVS 


1 - 




15(1 + ^73)' 


(/r(15 + p^)xo - 2(5 - M^/3)yo) + 


(109) 


The error estimation O {QL/ = O (1/SNR) comes from taking into account the next term in the T'd expansion of the 
complete action expression ( l 66 t . The estimation O {'y^QL^p^) appears as the estimation of the contribution of the nonlinear 
biquadratic terms, see Eq. (Unil, originating from the third line in the action (| 6 ^ . Indeed, it is obvious from the expression 
(l98b that any extra power of the field u{z) results in an extra factor y'QL. Substituting the leading order expression for fi'cZ 
in the third line of Eq. (l 66 b for these terms we arrive at the estimation O (yy^QL^p^y That is why the formal parameter of 
the perturbation theory for the nonlinear terms is QLj^L^p^. Note the quantity QLj^L^p^ = N^^Lp'P is the very parameter 
determining the upper limit of the intermediate power region (l23b . 

Finally, from Eq. (IHFb for the exponent factor and from Eq. ( 1109b for the pre-exponent factor we arrive at the expression 


exp 


P[Y\X] = 


(1 + 4^^/3)xg - 2pxoyQ + yl 
QL{l + p?/i) 
■nQL^Jl + pP/‘i 


1 - 


p/p 


15(1 + ^73)' 


(/r(15 + p^)xo - 2(5 - M^/3)yQ) - 


p/p 


135QL(1 + /r73)' 


3 (57 + 33/^2-h 30) J 


-|p (47 -f 157 + 225) teg + (237 + 255/r^ - 90) xgyo + /i (20/ + 117 p.^ - 45) XoVo “ 


O 


if) 


0(7“<3LV") 


Now it is easy to show, that the normalization condition 


J DYP[Y\X] = 1 


( 110 ) 


( 111 ) 


is fulfilled. 


C. Calculation of Pout\Y]. 

Let us consider the integral Pout[Y] = / VXPx\X]P\Y\X]. In our case the measure VX = dxdy, where x = Re{X}, 
y = Ira{X}, so we should consider the integral: 


dxdyPx[x,y]P\Y\X] . 


( 112 ) 


In the integral the scale of variation of the function Px[x,y] is P ^ QL. The scale of variation of the function P[y|2f] is 
QL, and this function has the form Eq. (l22b . therefore we can use Laplace’s method. To demonstrate that one can see that the 
function P[F|2f] depends on |2f|, xq = Re{X{Y e~''^ — X)/\XW, yo = Im{X{Ye~''^ — 2f)/|X|}, and reaches the maximal 
value at the point xq = yo = 0. Let us change the integration variables x, y to 771 , 772 , where rj = rji +7772 = {XP"^ — Y)e~^^'' '. 
Here is the phase of the Y. The inverse transformation reads X = {p + In the new variables 

the function P[U|X] reaches maximum at the point 771 = 772 = 0. The integral (II 12b takes the following form 

00 

J dpidp2Px i?e{(77+|r|)e-'^'^l''+l^ll'+''^‘'''},/777{(?7+|y|)e-'^-^l''+l^ll'+'‘^‘''’}] P[r|2s:], (113) 
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here we have used the fact that the Jacobian determinant for the variables transformation is equal to unity. Since P[y|X] 
reaches its maximum at the point p = 0 we can expand the functions Px[X] and P[F|X] in the vicinity of the point: 


(Px[-Re{Fe“*^},+ terms proportional to 77 + ...) , 


P[Y\X] 


TrQLy'l + fi^ 


exp 


(1 + 4/i^/3)7?^ - 2/27717?2 + vl 

QL(l + /i 2 / 3 ) 


1 + terms proportional to 77 and—— 

QL 


(114) 

(115) 


where we have used the fact that in the vicinity of the point rj = 0 we have xq = —rji and 7/0 = —V 2 up to higher powers of 
77 . In Eqs. (II141) and ( II151 ) we have the parameter fl = 'yL\Y\'^. 

One can see that at large p the exponent contains three different terms: 

(1 + 4 / 12 / 3 ) 77 ^ - 2 p 77 i 772 + 77 I _ 4^ _ g PiP 2 3772 - 9772 
QL{l + fP/3) QL QLfl QLp? 


(116) 


Therefore to use Laplace’s method we have to transform our quadratic form 

. . , .T (l + 4/iV3)?7?-2/l77i772 + 77| 

{m,V2)A{r]i,V2) = — 


to the canonical form. The matrix of quadratic form is: 

1 


A = 


QL(l + p2/3) 


Qi(l +7(2/3) 


, , 4/i2 _ 

1 + — -M 


-L 


1 


The eigenvalues of the matrix A are 


Ai — 


QL 


1 + /( 


A + \/9 + 4p2 
3 + /i2 


1 / p - y/9 + 47i2 

A2 = 7TT I ^ + L- 


QL 

One can see that Ai 2 > 0, and at large fl they have the form: 


3 + /i2 


A ^A A 

QL' ^ 4gL7i2 ■ 


(117) 

(118) 

(119) 

( 120 ) 

( 121 ) 


Therefore at large ji ~ fl there are two parameters in the Laplace integral, one parameter is 1/QL, the other is \/{QLfl’^). 
To use Laplace’s method for the integral Lq. (II121 ) we have to impose two conditions P ^ QL, and P ^ QLfl^. These 
conditions lead to the two dimensionless parameters for Laplace’s method 


SNR»1, (122) 

i-f^QL^P)-^ » 1. (123) 


To calculate the integral Lq. (II 13b in the leading order in the parameters 1/SNR and (j^QL^P) we substitute the first term 
of the expansion Lq. (II 14b and the first term in the brackets of the expansion Lq. (II 15b to the integral Lq. (II 12b . After 
straightforward calculation we obtain: 


Px 


00 

J dr]idr]2P[Y\X] Px 


(124) 


To calculate corrections to the integral in parameters 1/SNR and 'y^QL^P we should take terms which are proportional to 77 
and 77 ^ in the product of expansions Lqs. (II 14b and (II 15b . Lormally the first correction to the integral should be of order of 
l/VSNR and yJ"f‘^QL^P, but it is zero due to the symmetry 77 —77 (the exponent contains only even combination of 77 ). 

Therefore we can estimate the order of the first nonzero corrections as C>(1/SNR) and 0{'j‘^QL^P). Therefore the result for 
the integral Lq. (II 12b can be written as: 


00 



dxdyPx[x,y]P[Y\X] = Px 




(1 + 0 (1/SNR) + O {-i^QL^P)) . 


(125) 


— 00 
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